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Abstract 

We study the charmed meson pair (D^D^ and D*D~) production m pp annihilation within an effective 
Lagrangian model that has only the baryon-meson degrees of freedom and involves the physical hadron 
masses. The reaction amplitudes include terms corresponding to the t-channel Kf and baryon 
exchanges and the ^-channel excitation, propagation and decay of the 'P(3770) resonance into the charmed 
mesons. The initial- and final-state distortion effects have been accounted for by using a simple eikonal 
approximation-based procedure in the same way as was done in our previous study of the pp ^ 
reaction within a similar model. The production reaction is dominated by the baryon exchange 
process and the corresponding total cross sections are predicted to be in the range of 0.18-0.7 pb for an¬ 
tiproton beam momenta varying between threshold and 20 GeV/c. The T'(3770) resonance contributions 
have a large influence on the differential cross sections of the D~D^ production reaction. 

PACS numbers: 13.60.Le, 14.40.Lb, 11.10.Ef,13.85.Fb 
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I. INTRODUCTION 


The first discovery of a charm-anticharm (cc) bound state (7/<A) [llllj] was made more than 30 
years ago. Yet a substantial part of the charmonium spectrum is still to be precisely measured. 
Due to several reasons, charmonium states (and other heavy quarkonium states) have played an 
important role in our understanding of quantum chromodynamics (QCD), the fundamental theory 
of the strong interaction. Within the range of momentum exchange in bound cc systems, the value 
of the strong coupling constant as is not so large to invalidate the application of the perturbative 
methods. Thus these states provide a unique laboratory to explore the interplay between perturba¬ 
tive and nonperturbative effects in QCD. The relatively small binding energy of the charmonium 
as compared to the rest mass of its constituents allows its description by the nonrelativistic 
proaches that simplify and constrain the analysis of the nonperturbative effects (see, e.g. Ref. [j], 
for a recent review). These mostly analytical methods are of considerable help in making sig¬ 
nificant progress in lattice QCD calculations, which have become increasingly more capable of 
dealing quantitatively with the nonperturbative dynamics in all its aspects (see, e.g.. Refs. MlSfl). 

Therefore, there is considerable interest in investigations of the production of charmonium 
states. Experimentally, they have been studied mainly in electron-positron and proton-antiproton 
annihilation processes. However, there are distinct advantages in producing cc states in the latter 
method where all the three valence quarks in a proton annihilate with their corresponding antiquark 
partners in an antiproton. This does not set any constraint on the quantum numbers of the final 
states enabling one to reach all the charmonium states by the direct formation. On the other hand, 
in electron-positron annihilation, the direct creation of final charmonium states is constrained to 
the quantum numbers of the photon = 1““). Other states can be reached only indirectly by 
other mechanisms. 

The PANDA (’’antiproton annihilation at Darmstadt”) experiment will use the antiproton beam 
from the Facility for Antiproton and Ion Research (FAIR) colliding with an internal proton target 
and a general purpose spectrometer to carry out a rich program on the charmonium production in 
proton-antiproton annihilation. The entire energy region below and above the open charm thresh¬ 
old will be explored in these studies. Charmonium states above the open charm threshold will 
generally be identified by means of their decays to DD |^8], unless this is forbidden by some 
conservation rule. 

The reliable estimation of the rates of pp ^ D°D° and pp D~D^ reactions (to be together 
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referred to as the pp ^ DD reaetion) at the PANDA energies is required for the aeeurate deteetion 
of the eharmonium states above the DD threshold. In addition, it is also important for other 
studies sueh as open eharm speetroseopy, the seareh for eharmed hybrids deeaying to DD, the 
investigation of the rare deeays and of the eharge-eonjugation-parity violation in the D-meson 
seetor. All these topies are the major eomponents of the PANDA physies program [1^. The 
aeeurate knowledge of these reaetions is also the primary requirement for investigating the ereation 


1 ih and eharmed D-mesie 


of the exotie flavored nuelear systems like eharmed hypemuelei |19- 
nuelei 112-14], 

The eross seetions of the pp DD reaetion have been ealeulated by several authors employing 
a variety of models. In Ref. lilSll . a nonperturbative approaeh has been used, whieh is based 
on the 1/A expansion in QCD, Regge asymptoties for hadron amplitudes, and a strin g m odel. 


Similar types of models were used in the ealeulations reported in Refs. [ilbL IlTh. In Ref. lllSh . the 


pp ^ D^D^ reaetion has been deseribed within a double handbag model where the amplitude 
is ealeulated by eonvolutions of hard subproeess kernels (representing the transition uu ^ cc ) 
and the generalized parton distributions, whieh represent the soft nonperturbative physies. This 
approaeh was earlier used in Ref. [Il9n to deseribe the produetion of N~N^ in pp annihilation and 
it resembles the quark-diquark pieture that was employed by this group to make predietions for 
the eross seetions of the D D^ reaetion in Ref. OOll . 

Reeently the produetion of DD in antiproton-proton annihilation has been studied within the 


Julieh meson-exehange model in Ref. [l2lh . 


the pp AA 1122, 


23D and pp A/A/ [124 


lis approaeh was employed earlier to investigate 


2511 reaetions. In this model, these proeesses are 


eonsidered within a eoupled-ehannels framework, where the initial- and final-state interaetions 
are taken into aeeount in a rigorous way. The reaetions proeeed via the exehange of appropriate 
mesons between p and p leading to the final baryon-antibaryon states. 

Apart from Ref. lllSn . where the ealeulated total eross seetion (crtot) for the pp D^D^ reaetion 


was reported to be below 10 nb, in the majority of the ealeulations, the magnitudes of the cr,or for 
this reaetion lie in the range of 10 - 100 nb. However, the predietions of various models differ 
signifieantly for the eross seetion of the pp D~D^ reaetion. 

In this paper, we present the results of our investigations for the eross seetions of pp ^ D^D^, 
and ^ ^ D~D* reaetions within a single-ehannel effeetive Lagrangian model (see, e.g.. Refs. 
ll26l - l28n i. where these reaetions are deseribed as a sum of the t-ehannel A/, S/, 2/^ baryon ex¬ 
ehange diagrams [see. Figs. 1(a) and 1(b)] and the 5-ehannel exeitation, propagation and deeay into 
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FIG. 1. Graphical representation of the model used to describe the p+p ^ (a) and p+p —> D~+D^ 

(b) reactions via t-channel exchange of charmed baryons. In (a) AJ and Sj in the intermediate line represent 
the exchanges of and baryons, respectively while in (b) represents the exchange of E^"'' baryon. 


the DD channel of the 'P(3770) resonance (presented diagrammatically in Fig. 2). The f-channel 
part of the model is similar to that used in our previous calculation 11291] of the pp ^ reaction 
that proceeds via the f-channel and meson-exchange processes. 

In the next section, we present our formalism. The results and discussions of our work are 
given in Sec. III. Finally, the summary and the conclusions of this study are presented in Sec. IV. 


II. FORMALISM 


The exchanges of both A^ and baryons contribute to the f-channel amplitude of the p + p ^ 
reaction. However, the process p + p ^ D~ + is mediated only by the exchange 
of the baryon in the f-channel. On the other hand, the 5-channel excitation, propagation and 
subsequent decay of the intermediate resonance state i/r(3770) contribute to both these reactions. 

To evaluate amplitudes for the processes shown in Figs. 1 and 2, we have used the effective 
Lagrangians at the charm baryon-meson-nucleon vertices, which are taken from Refs. [I29l-l33n. 
For the vertices involved in the f-channel diagrams we have 
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FIG. 2. The Feynman diagram to describe the p+p —» D+D reaction via ^-channel excitation, propagation 
and decay of the 'P(3770) resonance. 


-C.NBD - iSNBD^Niy^^D^B + H.C., (1) 

where and are the nueleon (antinueleon) and eharmed baryon fields, respeetively, and is 
the D-meson field. gNBo in Eq. (1) represents the vertex eoupling eonstant. 

The amplitude of the diagrams given in Figs. 1(a) and 1(b) is given by 


A(5) = i- 
T 


Snbd _ 

- {niB - irB12y 


^p(.kp) (Tp/ + ’^b) y^ ikpikp). 


( 2 ) 


where B represents the exehanged charmed baryon. q, ms and Fb are the momentum, mass and 
the width of the exchanged charmed baryon, respectively. The term that contains these quantities 
comes from the propagator of these baryons. The widths of the charmed baryons are taken from 
the latest Particle Data Group estimates 11341] . The coupling constants gNBo are adopted from Refs. 


m, 


3l|], as gNA+D - 13.50, gATL+D — 2.69 and 2.69. From these values it is expected that 


will dominate the t-channel production amplitudes. 


The olf-s 


Refs. [126, 


lell behavior of the vertices is regulated by a monopole form factor (see, e.g., 


27l|) 


FiiqB.) 




m 


B, 


Af 




( 3 ) 


where qp, is the momentum of the zth exchanged baryon with mass mg., d, is the corresponding 
cutoff parameter, which governs the range of suppression of the contributions of high momenta 
carried out via the form factor. We chose a value of 3.0 GeV for d, at all the vertices. The same 
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Ai was a 


so used in the monopole form faetor employed in the study of the reaetion pp —> A 


+ 


in Ref. within a similar model. It may be mentioned here that in the Jiilieh meson-exehange 
model ealeulations of the reaetion p + p ^ D + D presented in Ref. Olll . a form faetor of the 


following type has been used: 


FiiqB.) = 




( 4 ) 


(/l;)4 + (qj - ml)^ 

with a d, of 3.5 GeV. This form faetor gives more weight to the lower momentum transfers. As 
diseussed below we use this type of form faetor at the resonanee produetion and deeay vertiees. 

In order to evaluate the diagram of Fig. 2, the effeetive Lagrangians are required at '¥pp and 
'TDD vertiees, whieh are written as 


= S'fpp l^P (yp + 


( 5 ) 


and 


-F'J>dD - S'¥DD (®D O/)) 


( 6 ) 


In Eq. (5) M represents the nueleon mass and 6^ is the 'T resonanee field. gxf>pp and are the 
eoupling eonstants at the 'Vpp vertex. Similarly, in Eq. ( 6 ) is the eoupling eonstant at the 
'TDD vertex and O 5 and represent the D and D eharmed meson fields. The values of the 
eoupling eonstants gvj/£,o£,o, g>PD-D+ and g^pp have been determined from the branehhm ratios for 
the deeay of 'T(3770) resonanee to the relevant ehannels as given in Refs. and |36]. We take 
g-j/^oflo = 17.90 (see also Ref. DTIl f. g>PD-z)+ = 14.10 and g^ipp = 5.12x10“^. The value of /cvy is 


fixed by fitting the eross seetions of the p + p ^ 'T(3770) ^ D°D° ealeulated within the effeetive 
Lagrangian model to that obtained within a semielassieal resonanee produetion and deeay model 
where experimental widths for the deeay proeesses 'T ^ pp and 'T ^ DD are used. This is 
diseussed in the next seetion. 

The amplitude of the proeess p + p ^ 'T(3770) ^ DD (Fig. 2) is written as 


A('I') = [Mr, + ^ (h - kof. (7) 

where kp and kp are the momenta assoeiated with the final-state D and D mesons, respeetively, 
and Sinv is the square of the invariant mass assoeiated with the 'T resonanee. It may be mentioned 
that the denominator of the 'T propagator leads to a eross seetion that has a pole in the vieinity 
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of the physical mass value (mvp) of the 'P resonance, which is taken to be 3773.15 MeV. The 
total width (Tvj/) of this resonance is 27.2 MeV Odil . A similar approach for the denominator 
of the 'P propagator has also been adopted in Refs. OSl - koll . This procedure is inspired by the 
extreme vector-meson dominance hypothesis where contributions of the vector-meson resonance 
are included in the vicinity of the relevant kinematical regime and in the far off-shell vector-meson 
kinematical regions it is considered as background that is generally quite weak. 

In order to account for the off-shell effects due to the internal structure of the intermediate 
charmonium states, we introduce vertex form factors for the 'P resonance. In our calculations the 
shapes of these form factors are given by Eq. (4) with cutoff parameters 


Tip - mip -I- qAqcd, 


( 8 ) 


where Aqcd = 240 MeV, and a is an adjustable parameter. We have used a = 7.5 for amplitudes 
involving the 'P(3770) resonance. 

A widely used approach is to parametrize the total cross section for the process pp 


ij/(3n0) DD in a Breit-Wigner form (see, e.g. Refs. 1135 . 


M 


m) 


O-n 


pp-^ifi(3nQ)^DD 


27>j/ -l- 1 47r Sinv^Wi3nQ)^pp rxj/(3770)-^Z)D 

(2jp + l)i2j, + 1)3 


-2 

tot 


( 9 ) 


ipp i^inv ‘^myT, 

where 7vp is the spin of the resonance and qpp is the momentum in pp channel. rij/( 3770 )^pp is 
the partial width for the production process pp ^ T* and rxj/( 3770 )^pD is the partial width for the 
'P ^ DD decay. Ytot is the total width of the 'P resonance. Attempting to account for the possible 
self-energy contributions to the formation and decay of the resonance, some authors (see, e.g .. 
Ref. [I 42 I 1 I have introduced an energy dependence to the width rij,( 3770 )^p£) adopted from Ref. ll43n . 
The ansatz for this energy dependence involves the range of DD interactions, which is treated as 
a free fitting parameter. For the sake of simplicity and in order to keep the number of adjustable 
parameters as small as possible, here we refrain from such a more elaborate approach. We use 
a constant width, which is a good approximation in view of the very narrow line width of the 
'P(3770) resonance. In fact, close to the resonance pole, the expression given by Eq. (9) is indeed 
a very good approximation of the exact result ll40ll . 

In our numerical calculations, we have used for Y^i( 3 iio)^dd the values extracted from the 
branchhm ratios of this decay as given in the latest compilation of the Particle Data Group (PDG) 
group [l34n . The width rxi/(3770) 


pp is obtained from the branching fraction = 7.D29XlO ^ 
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as reported in Ref. 11361] . We have taken the width eorresponding to the upper limit of whieh 


was also employed in the determination of the eou pling eonstant g^pp used in Eq. (7). 

From the studies of the A“A^ produetion 1124 1291]. it is well known that the magnitudes of 
the eross seetions depend very sensitively on the initial-state distortion effeets. In faet, the pp 
annihilation ehannel is almost as strong as the elastie seattering ehannel. This large depletion of 
the flux ean be aeeounted for by introducing absorptive potentials that are used in optical models 


or in coupled-channels approaches [122 - 




M 


451] ■ In this work, instead of employing such a 


detailed treatment, we use a procedure that was originated by Sopkovich [|46|] and was employed 
in Ref. for describing the pp A“A^ reaction. In this method, the transition amplitude for 
the reaction pp ^ DD with distortion effects is written as 


rj-pp^DD _ -^q^ppJ'^P^DD 


( 10 ) 


where is the transition matrix calculated within the plane-wave approximation and QpP and 

the are the operators describing the initial- and final-state elastic interactions, respectively. 

For the present purpose, we neglect the real part of the baryon-antibaryon interaction. Con¬ 
sidering the pp initial-state interaction (ISI), we describe the strong absorption by an imaginary 
potential of Gaussian shape with range parameter p and strength Vq. By using the eikonal ap¬ 
proximation, the corresponding attenuation integral can be evaluated in a closed form. Similar to 


Refs, m, 


47h . we obtain for QI’p 


QpP = exp{ — exp{-p^b^)\, 
pk 


( 11 ) 


where b is the impact parameter of the pp collision. E and k are the center-of-mass energy and the 
momentum of the particular channel, respectively. In our numerical calculations, we have used 
the same values for the parameters Vq, p and b as in Ref. It may be noted that with these 
parameters we were able to get cross sections for the pp ^ AA strangeness production reaction 
in close agreement with the corresponding experimental data. Furthermore, our cross sections 
for the pp A“A^ reaction were similar in magnitude to those reported in a coupled-channels 
meson-exchange model calculation [l24l] . 

In the case of the pp A“A^ reaction, it has been shown in Ref. [124 that, because of the 


strong absorption in the initial channel, the production cross sections were rather insensitive to 
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FIG. 3. The distortion factor Vlpp [defined by Eq. (11)] as a function of the antiproton beam momen¬ 
tum. The values of parameters //, Vq and b were taken to be 0.3369 GeV, 0.8965 GeV and 0.3270 GeV“\ 


respectively, which are the same those used in Ref. 12911 . 


the final-state interaetion (FSI) between and A^. In Ref. 112 lh . the effeets of DD FSI were 
investigated by a ppro ximately extending in the charmed meson sector the nn KK model of 


the Jiilich group I148h . However, these calculations have sizable uncertainties and even then the 


effect of FSI is not big. In our procedure, unlike the pp ISI, it is not possible to put any constraint, 
experimental or otherwise, on the choice of the DD FSI distortion parameters. Therefore, in order 
to keep the number of free parameters small, like our study on A“A^ production, we concentrate 
only on the initial-state interaction in this study. It should be mentioned that also in calculations 


reported in Refs. [115 


170 the meson-meson FSI effects were not considered. 


The distortion effects could lead to the reduction of the undistorted cross sections by several 
orders of magnitude depending upon the values of the parameters Vq, p and b. In Fig. 3, we have 
shown Q.pp as function of the beam momentum pp with the values of parameters p, Vq and b being 
0.3369 GeV, 0.8965 GeV and 0.3270 GeV“^ respectively. We see that Qpp increases gradually 
as Pp goes beyond the threshold and becomes almost constant at higher values of Pp - in this 
region the values of E and k are roughly equal. With this klpp, the undistorted cross sections would 
dampen by almost 2 orders of magnitude. As far as dependencies on the parameters are concerned, 
increasing Vq obviously increases the damping, while decreasing p has the same effect. 


9 













FIG. 4. Total cross section for the pp T'(3770) —> reaction as a function of the antiproton beam 

momentum. The solid line represents the cross sections calculated within the resonance production and 
decay model [Eq. (9)], while the dashed line shows the effective Lagrangian model cross sections [obtained 
with the amplitude given by Eq. (7)]. 

III. RESULTS AND DISCUSSIONS 


One can estimate the total cross sections for the pp and pp D~D^ reactions around 

the ‘F(3770) peak with the help of Eq. (9). With the values of various widths as extracted from the 
experimental data as specified in the last section, the total cross sections calculated by Eq. (9) can 
indeed be used to fix some of the parameters of the effective Eagrangian model (EEM). We deduce 
the parameter in Eq. 5 by comparing the total cross section for e.g., the pp ^(3770) ^ 
reaction calculated within the EEM [by using the amplitude given by Eq. (7)] with that obtained by 
Eq. (9). In the EEM calculations the pp ISI has been included. In Eig. 3, we show this comparison 
where the value of /cxp is taken to be 6.0. 

In Eig. 4, we see that the peak of the EEM cross sections coincides with that of the resonance 
model. Also, in the vicinity of the resonance peak the two models predict the same magnitudes of 
the cross sections and the same line shapes. The value of this cross section is about 12 nb at the 


peak position, which is close to the lower limit of a tot predicted in Ref. ll4flh . 


In Eigs. 5(a) and 5(b) we present the results for the total cross sections of pp ^ and 


pp ^ D D"" reactions, respectively, for antiproton beam momenta (p^^*) in the region of 6.4 - 6.8 
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FIG. 5. (a) Total cross section for the pp —> reaction as a function of the antiproton beam momentum. 

The solid and dashed lines represent the cross sections obtained by the coherent sum of the t-channel 
baryonic exchange and the s-channel 'P(3770) resonance excitation amplitudes, and by including the t- 
channel baryon exchange contributions only in the amplitude, respectively, (b) The same for the pp ^ 
D D^ reaction as a function of the antiproton beam momentum. The full and dashed lines have the same 
meaning as in (a). 


GeV/c, where the ¥(3770) resonanee is expected to show its impact. In these figures, solid lines 
represent the cross sections where the amplitudes of the f-channel baryon exchanges and the s- 
channel ¥(3770) production and decay are coherently added, while the dashed lines include only 
the f-channel baryon exchange contributions. We see in Fig. 5(a) that for the pp ^ reaction 
the cross sections corresponding to the f-channel baryon exchange processes are fairly large (~ 
300 nb) in the resonance peak (RP) region, so the inclusion of the ¥(3770) resonance contributions 
does not make any dramatic effect. It just produces a small kink in the cross section around the 
RP momentum. This is in contrast to the results of Ref. |^, where the ¥(3770) resonance led to 


an enhancement of almost a factor of two in the cross sections of the pp reaction around 
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the RP region. The reasons for this difference can be attributed to two facts. First, our t-channel 
baryon exchange cross sections are larger than those of Ref. ll40n (about 300 nb as compared to 
only about 40 nb) in the RP region, and second, our cross sections for the 'P(3770) resonance 
excitation are smaller than those of Ref. ll40ll in this region. 

On the other hand, in our study the effect of the ¥(3770) resonance is quite prominent for the 
pp D D^ reaction as can be seen in Fig. 5(b). In our model, the t-channel baryon exchange 
contributions to the cross sections of this reaction are strongly suppressed. The reason for this 
is that only the 2^^ exchange mediates the t-channel amplitude in this case, which becomes very 
small because of the much smaller coupling constant and somewhat larger mass of the exchanged 
baryon. The ratio of the absolute magnitudes of the pp and pp D~D^ reactions is 

roughly proportional to (gNAto/gN^ ++ d)"*, which leads to a reduction in the D D* production cross 
section over that of by nearly a factor of 650. 

It should, however, be mentioned here that in the consideration of the ISI within the distorted- 


wave Born approximation (DWBA) approach [Ref. [12111 1. two-step transitions of the form pp ^ 
hn D~D^, are generated. Because the AJ exchange can contribute to the hn D~D^ transition 
potential, this exchange is no longer absent. Therefore, these two-step mechanisms can enhance 
the D D^ production cross sections. Indeed, in Refs. Oil. koll the cross sections for D~D^ pro¬ 
duction are even larger than those of the production. On the other hand, such two-step 


mechanisms are out of the scope of our as well as of Regge model yj 


17h calculations. There¬ 


fore, in studies within these models the cross sections of the pp ^ D reaction are suppressed 
as compared to those of the pp ^ reaction. 

In Figs. 6(a) and 6(b) we show the differential cross sections (DCS) for pp and 

pp D D^ reactions, respectively, at the beam momentum of 6.57 GeV/c, which corresponds to 
the ¥ resonance invariant mass (= 3770.24 MeV. Therefore, these cross sections represent 
the angular distributions of the produced charmed mesons at practically the resonance peak. In 
Fig. 6(a), we note that the inclusion of the ¥(3770) resonance alters significantly the cross section 
obtained with only the t-channel baryon exchange term (dashed line). Strong interference effects 
of various amplitudes are evident. At the backward angles, the t-channel baryon exchange and 
i’-channel ¥(3770) amplitudes interfere destructively while at forward angles this interference is 
constructive leading to the strong forward peaking of the angular distribution. Because of this 
interference effect, the DCS do not exhibit the type of shape that is expected from a pure p-wave 
resonance dominated amplitude. 
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FIG. 6. (a) Differential cross section for the pp —> reaction at the antiproton beam momentum of 

6.57 GeV/c. The solid and dashed lines have the same meanings as those in Fig. 5. (b) The same as in (a) 
but for the pp —> D“D^ reaction. 

On the other hand, in Fig. 6(b), the differential eross seetion obtained by adding the 5-ehannel 
'F(3770) terms to the t-ehannel baryon exehange amplitudes shows a dominant p-wave type of 
angular distribution. This is due to the fact that contributions of the ^(3770) resonance term are 
significantly stronger than those of the t-channel baryon exchange term in this case. This was 
apparent already in Fig. 5(b). However, even though the t-channel baryon exchange amplitudes 
are relatively quite small, they do introduce some distortion to the angular distribution of the D 
mesons arising from the decay of the 'F(3770) resonance through the interference terms. This is 
evident from the asymmetry of the cos(0) distribution depicted by the solid curve. 

In general, the interference patterns in the differential cross sections depend quite sensitively 
on the relative magnitudes of the t-channel baryon exchange and 5-channel 'F(3770) amplitudes. 
Thus, they provide a critical check on the coupling constants that enter into these amplitudes. 
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FIG. 7. Total cross section for the pp reaction as a function of the antiproton beam momentum. 

The contributions of t-channel and baryon exchange processes are shown by dashed and dashed- 
dotted lines, respectively, while the dotted line represents the cross sections of the ^-channel T'(3770) res¬ 
onance excitation. The cross sections corresponding to the coherent sum of these amplitudes are shown by 
the full curve. 


Therefore, measurements of the differential eross seetions of these reaetions in future experiments 
would be useful in fixing these eoupling eonstants. 

Next, we diseuss the eharm meson produetion in antiproton-proton annihilation at higher beam 
energies. In Fig. 7, we show the total eross seetion of the p + p ^ reaetion for antiproton 
beam momentum varying in the range of threshold to 20 GeV/e. In this figure, the roles of various 
f-ehannel baryon exehange and the 5-ehannel 'T(3770) resonanee exeitation proeesses have been 
investigated. We note that cr,o, for this case is almost solely governed by the Ac-exchange mech¬ 
anism in the entire range of the antiproton beam momentum. The contributions of E^-exchange 
terms are lower by about 3 orders of magnitudes. This can be understood from the approximate 
proportionality of the ratio of atot of the two exchange processes to the fourth power of the ra¬ 
tio of the coupling constants of respective vertices involved in the corresponding amplitudes as 
discussed above. The ^-channel 'T(3770)-exchange term contributes negligibly to the CTtot of this 
reaction at higher beam momenta. 

We further note in Fig. 7 that the a,ot peaks at of about 9 GeV/c. This is in agreement with 


14 






lo"* 

110 -^ 

10 '"* 

10 '® 

6 8 10 12 14 16 18 20 

(GeV/c) 

FIG. 8. Total cross section for the pp —> D D^ reaction as a function of the antiproton beam momentum. 
The contributions of t-channel baryon exchange and ^'-channel 'P(3770) resonance excitation are shown 
by dashed and dashed-dotted lines respectively. The cross sections corresponding to the coherent sum of 
these terms in the total amplitudes are shown by the full curve. 



the results of the Regge trajeetory model ealeulations of Refs. [Il6ll and [Il5h . At the beam momen¬ 
tum of interest to the PANDA experiment (15 GeV/e), total eross seetion for the produetion 
reaetion predieted by our model is about 570 nb. This should be eompared with the results reported 
by other authors for this beam momentum. In Refs. [Q] and ll^ the eorresponding eross seetions 


are approximately 100 and 70 nb, respeetively, while in Ref. lllSh it is less than 10 nb. We reeall 
that in Refs. y and ll^ . the ISI effeets were ineluded by following an eikonal-model-based 
proeedure similar to that of our study. Of eourse, uneertainties in this method ean not be ruled out. 
However, in our ealeulations, we have taken the same ISI parameters in Eq. (II) as those used in 
our previous study ll29ll of the pp N~N* reaetion within a similar model. These parameters 
were eheeked by reprodueing the near-threshold eross seetion predieted within the Jiilieh meson- 


exehange model ealeulations of this reaetion reported in Refs. [l24l. 
treated more rigorously within a eoupled-ehannels method. 


2511 where ISI effeets have been 


In Fig. 8, we present the total eross seetions for the pp D~D^ reaetion as a funetion of the 
antiproton beam momentum. We see that the cr,o, of this reaetion is strongly suppressed eompared 
to that of Fig. 7. This was seen already in Fig. 5 at near-threshold beam momenta. A similar 
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FIG. 9. Differential cross section for the pp —> reaction at the antiproton beam momenta of 8, 10, 

12 and 16 GeV/c as indicated in the figure. The solid lines show the cross sections where the amplitudes of 
the r-channel A^, and baryon exchange and ^-channel T'(3770) resonance terms are coherently added. 
The dashed lines show the cross sections where the amplitudes include the contributions of the -exchange 
term only. 


suppression of D cross sections relative to those of has been noted in Refs 


- [Q], 


ttlTD, 


and llisll . This can be attributed to the much smaller S^^-exchange vertex coupling constant in 
comparison to that of the A^-exchange vertex. However, in the coupled-channels meson-exchange 
model, the initial-state inelastic interactions could enhance the D~D^ cross sections significantly 
as was discussed earlier. It would be quite interesting to test the predictions of various models for 
the cross sections in the PANDA experiment at the upcoming FAIR facility. 

As was already noted in Figs. 6(a) and 6(b), the differential cross sections provide more ex¬ 
plicit information about the reaction mechanism. These cross sections involve terms that weigh 
the interference terms of various components of the amplitude with the angles of the outgoing par¬ 
ticles. Therefore, in general the contributions of different mechanisms are highlighted in different 
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FIG. 10. (color online) Differential cross section for the pp D~D^ reaction for antiproton beam momenta 
of 8, 10, 12 and 16 GeV/c. The dashed lines show the contributions of the baryon exchange terms only 
while the solid line represents the cross sections where the amplitudes of baryon exchange and 'l'(3770) 
resonance terms are coherently added. 

angular regions. In Fig. 9, we show the predietions of our model for the differential eross seetions 
for the pp ^ reaetion at the beam momenta of 8, 10, 12 and 16 GeV/c. In this figure, the 
solid lines represent eross seetions that include the coherent sum of the amplitudes corresponding 
to the f-channel and baryon exchanges and the x-channel 'I'(3770) resonance terms, while 
the dashed lines show the cross sections where the amplitudes include contributions of the A^- 
exchange term only. Since, the 2^ baryon exchange and the 'I'(3770) resonance contributions are 
quite small compared to those of the A^ exchange, the interference effects of various terms in the 
amplitudes are not significant at higher beam momenta. We notice that with increasing beam mo¬ 
mentum cross sections are more and more forward peaked. This indicates the growing importance 
of the f-channel exchange terms with increasing beam momenta. 

On the other hand, in the differential cross sections of the pp ^ D~D^ reaction, the interference 
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effects of the t-channel baryon exchange and the 5-channel 'F(3770) resonance terms are 
visible even at higher beam momenta, as can be seen in Fig. 10. In this figure the dashed lines 
represent the cross sections when the contributions of only the baryon exchange terms are 
included in the amplitude, while the solid lines show the results where the amplitudes of 
baryon exchange and 'F(3770) resonance terms are coherently added in the cross sections. One 
notices that the inclusion of the 'F(3770) resonance terms changes the cross sections drastically 
at the backward angles. This effect is visible even at the beam momentum of 16 GeV/c. Thus 
the measurements of the angular distributions of the pp ^ D D^ reaction even at higher beam 
momenta can provide signals for the ¥(3770) resonance. Such a study would be complimentary 
to the methods proposed in Ref. [l49n 

It should be mentioned here that the physics of the charmed DD-meson production m pp an¬ 
nihilation for beam momenta in excess of 3 GeV may also involve vector resonances other than 
¥(3770). Some of these resonances are 7/¥, ¥(25), ¥(4040), and ¥(4160). However, only 
¥(3770) whose mass is just above the DD production threshold, has a substantial branching ratio 


(about 93%) for decay into the DD channel [1341] . Masses of both 7/¥ and ¥(25) are below the 


DD production thresholds and their widths are about 2 orders of magnitude smaller than that of 
¥(3770). Therefore, their decay to the DD channel is not possible even from the higher ends of 
their mass spectrum. Nevertheless, because the mass of ¥(25) (3.686 GeV) is just below the DD 
threshold, it may possibly decay to this channel due to off-shell effects. However, no branching 
ratio is known for this decay mode as per the latest PDG compilation 11341] . Therefore, we have not 
considered this resonance in our work. 

On the other hand, the masses of the resonances ¥(4040) and ¥(4160) are well above the 
DD threshold, and therefore they can decay to the DD channel. However, the branching ratios 
for these decays are hardly quotable according to the latest PDG compilation. Nevertheless, in 
Ref. [50] the branching ratios for the decays ¥(4040) ^ DD and ¥(4160) ^ DD have been 
estimated by fitting to the e^e~ DD data of the Belle Collaboration [15 ll] in the invariant mass 


region of 3.8 - 4.3 GeV. In this procedure these resonances are parametrized in terms of the Breit- 
Wigner form. The estimated branching ratios for the two decay channels are found to be (25.3 
+ 4.5)% and (2.8 + 1.8)%, respectively. These can be used to obtain the coupling constants 
5 t(4040)5d and gT( 4 i 60 ) 5 z)- However, to calculate the cross sections fox p + p ^ ¥(4040) ^ DD 
and p + p ^ ¥(4160) ^ DD processes, we also require the coupling constants gvp( 4040 )pp and 
gxj/( 4 i 60 )pp, about which no information is available. In any case, we performed calculations for 
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the cross sections of these reactions by taking the values of the gT( 4040 )pp and g'p{ 4 i 60 )pp to be the 
same as that of gT( 3770 )pp- We find that the resulting cross sections in the relevant region are much 
smaller in comparison to the total cross sections shown in Figs. 7 and 8. Therefore, the inclusion of 
resonances 'T(4040) and T'(4160) hardly lead to any noticeable change in the overall conclusions 
of this paper. 

Finally, we discuss the uncertainties and the range of the validity of our results. The theoret¬ 
ical approach (the effective Lagrangian model or ELM) considered in this work has mesons and 
baryons as effective degrees of freedom. This model will be valid in the energy range where con¬ 
sideration of explicit quark degrees of freedom is not required. Therefore, our model is certainly 
applicable in the range of beam momenta (from threshold to 20 GeV/c) considered in this work. 
However, the calculations performed within this model are sensitively dependent on the values of 
the coupling constants at various vertices involved in the t-channel and i’-channel diagrams, on 
the shape of the form factor and the value of the cutoff parameter involved therein, and on the 
parameters involved in the initial- and final-state interaction scattering matrices. The extents of 
uncertainty in our results due to all these issues are discussed in the following. 


We have taken t 


from Refs, yl 


52L 


le coupling constants (CCs) at the vertices involved in the t-channel diagrams 


5311 where they have been fixed by using the SU(4) symmetry arguments in the 


description of the exclusive charmed hadron production in the DN and DN scattering within a one- 
boson-exchange picture. The same coupling constants were used in the description of the charmed 


hadron production within the Julich meson-exc 


range model in Refs. [124 


40h . Furthermore, these 


coupling constants were also used in Ref. BOh to investigate the role of intrinsic charm in the 


nucleon using a phenomenological model formulated in terms of effective meson-baryon degrees 
of freedom. Thus, the coupling constants used in the calculations of the t-channel diagrams of 
our model are quite standard. The CCs at the vertices involved in the i’-channel diagrams are 
determined from the experimentally determined branching ratios of the decay of the ^^(3770) 
resonance into the relevant channels. Therefore, uncertainties in our cross sections due to the 
coupling constants are minimal. 

There may indeed be some uncertainty in our cross sections coming from the shape of the form 
factor [and the value of the cutoff parameter (T,) involved therein] that are used to regulate the 
off-shell behavior of various vertices. As stated above, we have employed a monopole form factor 
as given by Eq. 3, with a T, of 3.0 GeV. A form factor of a different shape and/or a different value 
of the cutoff parameter would lead to a different cross section. Eor example, using a quadrupole 
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form factor [of the type given by Eq. 4] with the same eutoff parameter leads to enhaneement in 
the eross seetions by faetors of 3 - 4. Changing 4, from 3 to 3.5 inereases the eross seetions by a 
faetor of up to 2. We have tried to minimize these uneertainties in the eross seetions by using the 
same shape (monopole) of the form faetor and the same value of d; that were used in our previous 
study of the eharmed baryon produetion [l29n . As in other eases, issues related to form faetor will 
be finally settled onee the data beeome available on the eharmed meson and baryon produetion in 
pp annihilation from the PANDA experiment. 

The initial- and final-state interaetions, whieh are the important ingredients of our model, 
provide another souree of uneertainty in our results. We treat these effeets within an eikonal- 
approximation-based phenomenologieal method. Generally, the parameters of this model are eon- 
strained by fitting to the experimental data. Beeause of the laek of any experimental information, 
it is not yet possible to test our model thoroughly. The absolute magnitude of our eross seetions 
may have some uneertainties due to this. Nevertheless, in our study we have used the same set of 
distortion parameters that were used in our previous ealeulations of the eharmed baryon produe¬ 
tion in the same reaetion. These parameters reproduee the data for the AA ehannel and the eross 
seetions for the A^Ac ehannel ealeulated within the Jiilieh meson-exehange model where distortion 
effeets are treated more rigorously within a eoupled-ehannels approaeh. Therefore, the initial and 
final ehannel distortion effeets ineluded in our model are eheeked against the other independent 
sourees. 


IV. SUMMARY AND CONCLUSIONS 


In summary, we studied the p + p ^ and p + p ^ D~D^ reaetions by using a single- 
ehannel effeetive Lagrangian model that involves the meson-baryon degrees of freedom. The 
dynamies of the produetion proeess has been deseribed by the t-ehannel A^, and baryon 
exehange diagrams and also the 5-ehannel exeitation, propagation and deeay of the 'P(3770) res- 
onanee. The initial- and final-state interaetions have been aeeounted for by an eikonal type of 
phenomenologieal model. The eoupling eonstants at the baryon exehange vertiees were taken 


from Refs, yi 


52L 


3l, 


whieh were the same as those used in the study of the p + p ^ 
reaetion at the similar vertiees in Ref. 11291] . The CCs at '^pp, and vertiees have 


been determined from the branehing ratios for the deeay of 'T(3770) resonanee into the relevant 
ehannels as given in Refs. 11351] and 11361] . The off-shell eorreetions at various vertiees have been 
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accounted for by introducing monopole form faetors with a eutoff parameter of 3.0 GeV. The 
same form faetor with the same value of the eutoff parameter was also used in our study of the 
p + p —> A“A^ reaetion [l29n . The parameters involved in the initial-state interaetion seattering 
matriees were also taken to be the same as those used in Ref. 

Sinee the eross seetions of the p + p ^ D~D^ reaetion are strongly suppressed due the smaller 
eoupling eonstants of the vertiees involving baryon exehange, the inelusion of the 'T(3770) 
resonanee produees a sizable enhaneement in the p + p ^ D~D^ eross seetions around the reso- 
nanee energy. However, their effeet is not so strong in ease of the p + p ^ reaetion where 
the baryon exehange eross seetions are quite large - they vary between 100-400 nb for antiproton 
beam momenta between 6.4 - 6.8 GeV/e. Therefore, the inelusion of the 'T(3770) resonanee in 
this ease produees only a small kink in the total eross seetion near the resonanee energy. 

On the other hand, the differential eross seetions for both the reaetions are affeeted in a major 
way by the 'T(3770) resonanee eontributions for antiproton beam momentum near the resonanee 
peak. In ease of the p + p ^ reaetion, the T* resonanee eontributions introduee sizable 

reduetion (enhaneement) in the DCS at the baekward (forward) angles. For the p + p ^ D~D* 
reaetion, the shape of the DCS ehanges drastieally by the inelusion of the T* resonanee term - it 
ehanges to a p-wave type of distribution from a s-wawe shape. This drastie shape ehange of the 
DCS ean perhaps be exploited in a dedieated experiment at the PANDA faeility to pin down the 
*^(3770) resonanee. 

At higher antiproton momenta, the total eross seetion of the p+p ^ reaetion is dominated 

by the eontributions of the A^ baryon exehange. The eross seetion peaks around Pp^ of 9 GeV/e. 
At a Pp^ of 15 GeV/e, whieh is of interest to the PANDA experiment, the total eross seetion of 
this reaetion is about 550 nb whieh is at least 5 times larger than the largest value of this eross 
seetion reported previously. Of eourse, previous ealeulations have used different types of models 
that invoke explieitly the quark degrees of freedom in their ealeulations, whieh may make them 
more adequate for energies higher than those of the PANDA experiment. Therefore, it is not trivial 
to understand the reasons for the large differenee seen between their eross seetions and ours. The 
future PANDA experiments at FAIR are expeeted to elarify the situation. 

Within our model the eross seetions of the p + p ^ D D^ reaetion are strongly suppressed as 
eompared to those of the p + p ^ D^D^ reaetion. This is due to the faet that the latter is dominated 
by the A^ baryon exehange meehanism, while the former gets a eontribution only from the 
exehange whose eouplings are mueh lower than those of the A^ exehange vertiees. However, in 
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the coupled-channels meson-exchange model of Ref. II21I1 . the p + p ^ D~D^ cross sections are 
even larger than the p + p ^ ones. This is a result of the coupled-channels treatment in the 
incident channel which accounts effectively for two-step inelastic processes involving ’’baryon 
exchange.” 

The differential cross sections of the p + p ^ reaction at higher values of are 

strongly forward peaked and are so strongly dominated by the contributions of the baryon 
exchange terms that the interference terms of various other contributions baryon exchange 
and resonance) become insignificant. However, in the case of the p + p ^ D~D* reaction, 
differential cross sections have significant contributions from the interference terms of baryon 
exchange and the'T resonance process even at higher antiproton beam momenta. In view of these 
results, it should be possible to pin down the 'P(3770) resonance contributions in these reactions 
in dedicated experiments in the relevant antiproton energy regions. 
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